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SUMMARY 

The linearized drag integral for bodies of revolution at super- 
sonic speeds is presented in a double-integral form which is not 
based on slender-body approximations but which reduces to the 
usual slender-body expression in the proper limit. With the aid 
of a suitably chosen auxiliary condition, the minimum-external- 
wave-drag problem is solved for a transition section connecting 
two semi-infinite cylinders. The projectile tip is a special case 
and is compared with the Von Kdrmdn projectile tip. Calcula- 
tions are presented which indicate that the method of analysis 
gives good first-order results in the moderate supersonic range. 

INTRODUCTION 

In making the slender-body approximation to the linear- 
ized supersonic-flow theory for bodies of revolution, a basic 
approximation leads to replacing the axial sowce distribu- 
tion with the cross-sectional-area derivative. Slender-body 
theory, therefore, becomes linear in the superposition sense 
for cross-sectional areas as well as for sources or fields in 
contrast to linear supersonic-flow theory which is linear in 
the superposition sense for sources or fields but not for areas. 
Maldng the slender-body approximation, however, eliminates 
a large part of the Mach number dependence of the results. 
Lighthill (ref. 1) has shown that this basic approximation, 
for sufficiently smooth bodies, has the same mathematical 
order of accuracy as the linearized supersonic-flow equation. 
Ward (ref. 2) has extended the generality of slender-body 
theory by presenting a drag expression which is valid for a 
body with a finite slope at the base. Lighthill (ref. 3) has, 
at the price of a large increase in complexity, modified 
slender-body theory to include area-derivative discontinu- 
ities at a finite number of points. 

In 1935 Von K6rmdn (ref. 4) determined the minimum- 
wave-drag projectile tip. Later Sears (ref. 5) and Haack 
(ref. 6) determined minimum-wave-drag shapes for projectile 
tips and closed bodies of revolution subject to various com- 
binations of auxiliary conditions of constant length, constant 
caliber, and constant volume. The Von Kdrmdn tip and 
the Sears-Haack bodies are based on slender-body theory. 
In contrast to the minimum-drag bodies of given volume, 
which are accepted as reasonably good first-order results, the 
minimum-drag projectile tips and open-nosed bodies which 

i Supersedes NAOA TN 3189, 19M. 


have been derived from slender-body theory are subject to 
more severe restrictions on the range and conditions for 
validity. Busemann (ref. 7) pointed out this fact for the 
Von Kdrmdn projectile tip in 1941. In the late 1930’s 
Ferrari made studies of the minimum-drag projectile tip 
and length-caliber body problems. Recently, Ferrari (refs. 
8 and 9) has considered the length-caliber body and ducted- 
body minimum-drag problems on the basis of linear theory 
without resorting to slender-body approximations. He ob- 
tains source-distribution functions which involve elliptic in- 
tegrals that have to be treated numerically. 

This report presents the linearized-drag integral for bodies 
of revolution in a double-integral form which is not based on 
slender-body approximations but which reduces to Ward’s 
drag expression in the proper limit. With the aid of a suit- 
ably chosen auxiliary condition, the minimum-external-wave- 
drag problem for a transition section connecting two semi- 
infinite cylinders is solved. The source distribution, the 
minimum-wave drag, the slopes at the ends of the section, 
and the radius at one intermediate point are obtained in 
terms of elementary functions. The entire shape is given 
in an integral form amenable to numerical evaluation. 

SYMBOLS 


c 

constant factor in source-distribution func- 
tion, ft/sec 

Gd 

drag coefficient 

V 

drag, lb 

f 

source-distribution function 

f 

derivative of f with respect to its argument 

M 

free-stream Mach number 

r 

radius in cylindrical coordinates 

R 

general radius of point on body 

Ri 

radius of downstream cylinder 

Ri 

radius of upstream cylinder 

A'(x), A"(x) 

first and second derivatives of cross-sectional 
area 

U 

free-stream velocity, ft/sec 

V r 

radial velocity, ft/sec 

v x 

axial velocity, ft/sec 

X 

axial distance in cylindrical coordinates 

So 

position of first disturbance on control surface 
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axial distance for general point on body 


Sf 


variation in derivative of source-distribution 
function 

p free-stream density, slugs/cu ft 

X arbitrary constant 

£,ij dummy variables of integration 

£o, ijo position on axis where sources begin 

4 > disturbance potential 

All distances are made dimensionless by measuring them 
in the units of length of the transition section. 

ANALYSIS 


DRAG INTEGRAL 


The linearized-flow equation for the case of axial symmetry 
is (ref. 4) 


b z 4 > 1 b<j> Ri 2>V_, 
dr^dr P bJ 


(1) 


where <f> is the disturbance potential, p 3 =M*— 1, and M is 
the free-stream Mach number. The general solution of 
equation (1) is (ref. 10) 

JV-^r cosh „)*=-£'' (2) 

The corresponding disturbance velocities are 
J7(x-»- cosh 

(3) 

and 

~^~v r {x,r)=~ J f'{x—fir cosh u)p cosh u du - 

_1 f *-* f®(x-Qdt fA 

rJ— rJ(x-Q 3 - 0 b^ K 


The general solution is interpreted physically as the dis- 
turbance potential of a distribution of axial sources. The 
source-distribution function /(£) is fixed by satisfying the 
boundary condition that the body surface be a stream 
surface. 


Applying the momentum and continuity theorems to the 
fluid within a cylindrical surface of radius Ri and with the 
symmetry axis parallel to the stream (hereinafter referred to 
as the control surface) and using the fact that free-stream 
conditions prevail at the two ends of the cylinder, x— <» and 
2= — oo t yield the usual expression for the drag (ref. 4) 

D=— 2 v J pRiVt{x,Ii^v r {x,Ii^dx (6) 

If a position x=L on the control surface is chosen so that, 
for x^L, either v T or v x vanishes on the control surface, and 
if a® is the position of the first disturbance on the control 
surface, equation (5) reduces to 

D=— 2 ?r ( pRiv x (x,Ri)v r (x,R{)dx (6) 


Substituting expressions (3) and (4) for the velocities into 
equation (6), where £o=vo—Xo—PRi, the position on the 
axis where the sources begin, yields 


T) o C L [ f T ~^ /'(Qfr-e# T f(y)dr, 1 

2 J-oLJfo V(*"0 , -wJU’>o ^(x-vf-pR?\ d 

( 7 ) 


or 


~ ~ fT rx ~ pRl ntrmx-mdn i, ? „ 

[_j£o J»o ^(x—Q 3 —fi 3 Ri 3 V( x—itf—fPRfj 

( 8 ) 


The basic slender-body approximation reduces the integral 
over £ in equation (7) to f(x). The use of equation (4) for 

• • di* t? • 

small values of r and the boundary condition ^ gives 


f(x) proportional to the axial derivative of the cross-sectional 
area. The slender-body approximation thus reduces the 
integral of equation (8) to a double integral (over x and 17). 
In principle, the integral in equation (8) may be reduced to 
a double integral over £ and r, by interchanging the order of 
integration and performing the integration over 2. When 
this process is attempted in equation (8), however, the 
integration over 2 involves elliptic integrals and no essential 
advantage is gained. 


However, the drag integral may be reduced to a double integral in the following manner. Equation (8) is rewritten 
with the dummy variables £ and >7 interchanged. Addition of the interchanged and original expressions and division by L 
yield 

t-> r . d v m 

A, J{x-&-pRMx-rif- 0 *R l * 

Interchanging the order of integration yields 

CL-/3S, pi 

mr(v)<% dr, 

( 0 t/v 0 i 


(x—£+x—y)dx 


MR: VFFWVF vY-m 
irf£-Ri 


(10; 
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where the lower limit of the integral over x is £+/Si?i or i}+ pR x , whichever is the larger. The results of the integration 
over x is (see the appendix; for the details of the integration) 


F (x-S+x-ri)dx r u . x 

(x-Q(x- v )-pRY 

=cosh -1 

-lf+0*. 

7+0*1 

(L-QiL-ti-fRY 

h+fiR^x-tf-PRY^x-vY-PRi 2 L 

7+0*1 

pRi(Z-v) 

pRi(£— v) 


and the drag becomes 


rL-pRl C'L-pRi 

d=*p\ mm cosh - 1 


(£-£)(£- q)-/W 
pRi(£— v) 


d£ dr) 


( 11 ) 


( 12 ) 


Equation (12) gives the drag of an axial distribution of sources/® which begins at £> and is subject to the condition that 
there is no momentum flow through the control cylinder of radius R x for x ^ L. 


In the slender-body approximation, equation (12) reduces 
to Ward’s drag expression as follows. Use of the expression 
for the radial velocity (eq. (4)) for small values of r and 

the boundary condition give 

j(x)=^-A'(x) (13) 


where A'(x) is the axial derivative of the cross-sectional area. 
For small values of R x , pR x may be neglected in the limits of 
the integrals in equation (12) and the cosh -1 factor may be 
approximated by 


cosh -1 


{L-QjLm-W 


PRi{Z—v) 


=log [(L~a(L~ v ))- 
log |S-i7|-log— (14) 


Simple manipulation, based on the assumption that A'(£) 
is continuous and A' (£ 0 ) = 0, gives the Ward drag formula 

f*f* A"(QA"(v) log d&v+ 


cylinders are of radius Ra and R X) respectively. The length 
of the transition section is taken as the unit in which all 
distances are measured. The cylinder of radius R x is chosen 
as the control surface. 

Figure 1 is drawn for the case R 3 <^R X . However, the 
statement of the problem is the same as for RA>R\ since only 
the appropriate changes in sign occur. Equation (12) for 
the drag becomes 


D 


J 1-0*1 r 1-0*1 

-0*, J-0*, 


mm cosh -1 




pRife— v) 


d£di i 
(16) 


The auxiliary condition for the minimum-wave-drag 
problem is the requirement that the points ®=1, r—R x and 
x=0, r—Rz be on the same stream surface. In terms of 
mass flow this condition requires that the free-stream mass 
flow between the control cylinder of radius R x and the up- 
stream cylinder of radius R* equals the mass flow out of the 
control cylinder between x—p(R x —Rz) and x=l; that is, 

pU-rriR-f—Rz)— f 2-KpR x v t (x^R i)dx— Cons tan t (17) 

or 


^ A'{L) JJ A"(£) log {L-m- 1 — IA\L)Y log & 

(15) 

MINIMUM-DRAG PROBLEM 

The problem of determining the shape of a transition sec- 
tion connecting two semi-infinite cylinders for minimum ex- 
ternal wave drag may be formulated in the following way. 
(See fig. 1.) The upstream and downstream semi-infinite 



P f 

j0(s I -B J ) J - 


; -p«. m(x-m 

0*2 V(x — £) 2 — FRY 


dx=^r {R Y~R Y ) = Constant 


(18) 


Equation (18) is a particularly convenient form of the 
auxiliary condition in that the source-distribution function 
is the only unknown involved. A direct integration of the 
usual boundary condition would lead to an auxiliary condi- 
tion involving, in addition, the unknown shape of the section. 
Without further approximations, the minimum-drag prob- 
lem would be unmanageable. It appears impossible to 
obtain an auxiliary condition of constant volume involving 
the source-distribution function as the only unknown. 

Interchanging the order of integration in the left member 
of equation (18) and integrating over x yield 

P"^ 1 /'©/(l-f) 2 -^! 2 d^=~ (RY-Rfi^C onstant (19) 

J-0*2 Z 
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The problem of minimizing the drag given by equation (16), subject to the auxiliary condition given by equation (19), 
can be expressed in the usual manner of the calculus of variations. The expression 


<?(f , Ei, R^=^^mn-n) cosh- 


( 1 - 


-8(l-?)-/3W 

pRiii—y) 


d&v+2\ f 1 'V(0VU 

J -pH* 


(20) 


(where constant factors have been absorbed into the arbitrary constant X) must have a zero variation corresponding to an 
arbitrary small variation 5/' in/'; that is, 


<?(/'+5/', R u £*)-<?(/', R u RJ= 0 

Combining two terms, after interchanging dummy variables and neglecting the term involving (6/') s , yields 

l(l-Q(l- g )-W 






dr 1 +'K^(l-Q i -0 1 R 1 2 d£—0 


( 21 ) 


Since the variation 5/' is arbitrary, 


J l-flEl 

-SRa 


f(y) cosh" 1 


q-ea-n)-W 


dij+ 


pRi(t—y) 

W(l-f) s -W=0 (22) 

An integration by parts, with / (— fiR 2 ) = 0, yields 


X 


1 -0R1 


J(y)dv_ 


es* (f -^-y/a-v) 3 -^ 


=X= Constant 


(23) 


The integral equation (23) is the same equation that appears 
in the minimum-induced-drag problem in lifting-line airfoil 
theory. 

The solution of the integral equation 


X 


' F(y)dy _ 




= Constant 


is 


F(y) 


Oy+q 


-J(b—y)(y—a) 
Therefore, the solution of equation (23) is 




(24) 


The condition that /(— pRi)=0 requires that C\ — (7/3 7? 2 
therefore, 

m=C\l(r,+W 2 ) (1+pRx-r,) (26) 

The source-distribution function for minimum wave drag is 
part of an ellipse. The constant O is fixed by satisfying the 
auxiliary condition, equation (19). The final result for the 
source-distribution function for minimum wave drag is 


/(v)= 


4U{R?-R?) 


ir(l — j822j-f-/3.Ra) (l-\-3fiRi~\- ft Ri) 


V( 7 )+/3A a )(l-(- f}R i — 17 ) 
(26) 


Figure 2, which has an arbitrary vertical scale, shows the 
source distribution corresponding to the transition section 
of minimum wave drag. 



Figure 2. — Source distribution (arbitrary vertical scale) corresponding 
to transition section of minimum drag. Ellipso of which the distri- 
bution is a part is shown by dashed curve. 


MINIMUM-DEAG SHAPE 


The expression for the minimum-drag shape is analogous to the auxiliary condition given by equation (17) except that a 
control surface of arbitrary radius R is used, where R is the radius of the section at x=X (see fig. 1) ; the shape is given by 


fUv(R 3 — Ri)— f 2 tt pRvr(x,R)dx 

JPiB-Rf 

= 2 „f x rr* tMt=ML~\d x 

JPlB-Bj V(Z— £) !— / 3 2 F 2 J 


(27) 


Interchanging the order of integration, integrating over x, and substituting the source-distribution function yield 


(3 3 (R 2 -Rfi= 


8/3 W-ig* 2 ) 


v(l— /S-Ri+fi-Ra) (1+3/32?!-}- /3J?j) J-pR^ 


X 


V(X _Q 3 _^ 


V(H‘fi-K 3 ) (l+j3l?i— £) 


on 


(28) 


Since all the radii in equation (28) are multiplied by /3, one calculation gives a family of shapes for which, at a given value 
of X, f}R(X) ^Constant. Thus the minimum-drag shape obeys the Gothert similarity rule (retf. 11). 
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At the general point (X, R), the integral in equation (28) involves elliptic integrals and cannot be evaluated in elementary 
form. However, at one particular value of X, corresponding to the center of the source-distribution ellipse, the integral 
reduces to elementary functions. Thus, for 

xr_l+0 Ri-PR* 

X g 

^ (pRS+pRfia+fiBi+Wtf-SpBSR,* 
p 2[(i+fm 1 +miY-W*Ri 2 +PB/)\ 

Equation (28) may be solved by numerical-graphical methods. By a change of variable of the form 

y l+d-Ri — PRi i l+ABi+fJRa 

2 h 2 C0S 9 


equation (28) is converted into the following form suitable for numerical integration: 


m , B « = jsaaiaj 


— j 8Ri l-h/tBi-h/iBs 


o)-pR*dB 

'^■»0) 


The slope of the transition section at the general point (X, R) is found by differentiating equation (28). The result is 


dm) P O C X ~P R 


— PRi 


f 

~RUj- 




V(H-/9A 2 )(1 +pR l ~Q ^{X-tf-pR' 


a +pR 1 -pR i 


! I v 2 

U j-pr j YG+^ft)(l+ABr-^ 


-f) d£ 


For the general point, equation (31) is not integrable in elementary form. However, at the downstream end of the section 
(X=l, j£=Ri) it may be integrated directly. The result is 


V d(fiR) l V JWi (1 +pR l -pR 3 )(l +PRi 

L dX Jx-i L V2 (/PRY-PRY) Vl +PRi 


-bPRyj-ZPRy) 

-{-pRi-\-pRt 


At the upstream end of the section, equation (31) may be evaluated by a suitable limiting process. For example, if R=Rt 
and X—e (where e is a small quantity), the integrals in equation (31) may be evaluated if e is neglected in comparison with 
quantities of the order of unity. After the evaluation, « is allowed to vanish, and the result is 


vdtom r 

L dX Jx-o [ ‘ 


-\lpRj (1 — pRi -j- PR ?) (1 -j- 3pRi~\-t3R%) 

V^Vl ~\~pRi~\~PRi (fYRY fYRY) _ 


For Ri—0 (the projectile tip), the front slope is equal to the slope of the Mach cone. The back slope is different from zero 
for 


MINIMUM WAVE DRAG 

The expression for the minimum wave drag is obtained by 
inserting the source-distribution function (eq. (26)) into the 
drag integral (eq. (16)). The expression for the drag of the 
transition section is 

2*-pUW--R 3 y . . 

u a-pRi+PRJQ-+WRi+f}Rd J 

and the drag coefficient (based on the area of the section 
projected on a plane normal to its axis) is 

a 4 ] RY-RY\ (3 « 

Ud (i-pRi+PRJQ-+WRi+PRi) v ; 

REVERSIBILITY 

The results for the reverse-flow case may be obtained, 
formally from the preceding direct-flow case by interchang- 


ing Ri and R?. If, in addition, X is replaced by 1— X, the 
X-station will be measured from the same end of the con- 
figuration. The transition-section-shape expression (eq. 
(28)) can be shown to be an invariant relation between X 
and R under the transformation 

Ri—>Ri 

Rf—^R-i 

X-*l-X. 

Since the present analysis is a linear treatment, the reversi- 
bility of the minimum-drag shape c ann ot be considered as 
established to an order higher than the first. 

The drag coefficient (eq. (35)) is not symmetrical in Ri 
and Aj. However, f?C D is symmetrical to the first order in 
the small quantities flRi and pR* or, if fiRi and p R? are not 



116 


REPORT 1213 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


small, p?C D is symmetrical to the first order in the quantities 
pRi—pRi and PR 2 —PR 1 . The minimum wave drag thus is 
the same to the first order in the appropriate small quantities 
for the direct-flow and reverse-flow cases. 

DISCUSSION OF RESULTS 

The drag integral (eq. (12)) presented in this report de- 
pends upon the source-distribution function in the same way 
(except for slightly different weights of the kernel functions) 
that the slender-body drag expression (eq. (15)) depends 
upon the cross-sectional-area derivative. Since bodies with 
area-derivative discontinuities can be generated by continu- 
ous source distributions, equation (12) gives finite drag for 
certain types of bodies (for example, bodies with shoulders 
forward of the base) for which the usual slender-body drag 
expression diverges. The drag integral presented herein is 
thus applicable to a larger class of bodies than the Ward 
drag expression. 

Some indication of the range of validity of equation (12) 
is given by a comparison of the values of drag calculated by 
various methods for several cones. Figures 3 (a), 3 (b), 
and 3 (c) present cone drag coefficients plotted against pR t 
for semivertex cones of 5°, 10°, and 15°, respectively. The 
drags are calculated by the following methods: 

Method A — The drag integral (eq. (12)) with mass-flow 
continuity (eq. (17) with 7^=0) used as the boundary 
condition. Since the continuity condition is used in deriving 
the drag integral, this method uses a consistent boundary 
condition. 

Method B — The solution of the linearized supersonic-flow 
equation for a cone (/' (£) = Constant) used with the exact 

boundary condition and the exact pressure 

equation. This method may be called the exact linear 
method. Van Dyke (ref. 12) has presented cone drags by 
this method. 

Method C — The exact characteristics method for conical 
flow as given by reference 13. 

Method D — The drag integral (eq. (12)) with the boundary 
dv v 

condition This method was used only for the 5° 

cone. Each of the figures also shows the drag of the mini- 
mum-drag projectile tip of the same thickness ratio as was 
determined in the present analysis (curve I in fig. 3). 

Figures 3 (a), 3 (b), and 3 (c) show that the drag integral 
with the continuity boundary condition (method A) gives 
cone drags closer to the exact value (method C) than the 
exact linear method (method B) for values of P R x less than 
approximately 0.30. From a comparison of methods A and 
D for the 5° cone, the divergence of the cone drag calculated 
by method A as pRi-*l.§ is seen to be caused by the use of 
the mass-flow continuity boundary condition. These results 
indicate that equation (12) with equation (17) as the bound- 
ary condition gives good first-order results for Mach numbers 
above" the transonic range and less than the Mach number 


at which the tangent of the Mach angle is less than approxi- 
mately three times the slope of the body, that is, a moderate 
supersonic range. 

Projectile tip (Ri=0). — Figure 4 presents two families 
of minimum-drag projectile-tip shapes (j3jK 1 =0.2 and 0.4) 
calculated by equation (30). The Von Kfirmdn projectile 



Method A 

Drag integral with mass-flow-continuity bound- 


ary condition 

Method B 

Exact linear method 

Method C 

Exact characteristics method from reference 13 

Method D 

Drag integral with 

I 

Drag coefficient of minimum-drag projeotllo tip 


(a) 5° cone. 

(b) 10° cone. 

(c) 15° cone. 

Figtjee 3. — Comparison of drag coefficients for three cones determined 
by several methods and for minimum-drag projeotllo tip of same 
thickness ratio as the cone. 
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tip and the limi ting case /3 jR,= 1.0 are shown for comparison. 
The Mach number and thickness-ratio dependence are 
related by the fact that the minimum-drag shape obeys the 
Gfithert similarity rule (ref. 11). The slope at the base 
varies from zero to the cone slope as /S.Ri varies from 0 to 1.0. 

In both the Von Kdrmfin analysis and the present analysis 
the basic assumptions of linearized theory are strongly 
violated at the nose of the projectile tip. Specific results of 
either analysis at the nose therefore cannot be accepted with 
confidence. However, inasmuch as linear theory often 
qualitatively predicts correct trends even outside its range 
of quantitative validity, it is interesting to note that in both 
analyses the nose of the minimum-drag projectile tip is as 
blunt as the analysis will allow (infinite slope for the Von 
Karmun tip and tangent to the Mach cone in the present 
analysis). It is probably true that the exact shape of a 
projectile tip for minimum wave drag has a blunt nose. 

Figure 5 presents the ratio of the minimum-wave-drag 
coefficient based on the base area to the Von Kdrmdn drag 
coefficient as a function of pR x . The ratio has a minimum 
value at pRi= 0.33 and diverges as pRi~*1.0. The results 



Figure 4. — Minimum-drag projectile-tip shapes. Von K&rmdn 
projectile tip is shown for comparison. 



Figure 6. — Ratio of minimum- wave-drag coefficient to Von KArmdn 
drag coefficient as a function of 


for the cone discussed previously suggest that the drag 
diverges as /JjRi— >1.0 because of the use of mass-flow con- 
tinuity (eq. (17)) for the boundary condition and that the 
results are reasonably good for pRi less than about 0.3. 

General case. — Figure 6 presents the radius at an inter- 
mediate point and the slopes at the two ends of the mini- 
mum-extemal-wave-drag transition section for four configu- 
rations with pR 1 —pR 2 =0.4 : . Figure 7 presents these three 
quantities for five configurations for which pR l —pR 2 =0.05. 
Equations (29), (32), and (33) indicate that for PRi— pR 2 <$ 

j „D i o-n ^(/3Bi— jSfij)(l-f/SjBi— jfitB*) , 

and pR,+pR 2 > i-^3(fiR 1 —pR ) shape is a re- 

verse curve lying within the truncated cone connecting the 
two cylinders for a distance downstream of the nose which 
increases as pRy—pR 2 decreases. The entire shape was 
calculated by equation (30) for the case — jS7? 3 =0.05 and 

pRi+f}R i = 0.118, which corresponds to the minimum front 
slope for pR x — pR a — 0.05 (shown in fig. 7). This reverse- 
curve effect is more pronounced in the range of small values 
of pRi—pRs, the range where the linearized treatment is 
certainly correct. 

For a comparison with the general case, the drag of 'a 
truncated-cone section for pRi—pRi=0.05, pRi+pRn=0.118, 
and P—1 was calculated by the method of reference 14. 
The truncated-cone drag coefficient (0.0228) is 18 percent 
higher than the minimum drag coefficient (0.0193). 



Distance along axis, X 


Figure 6. — Radius at intermediate point and slopes at two ends of 
minimum-external-drag transition section for four configurations. 
0.4. 
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CONCLUDING REMARKS 

The drag integral for bodies of revolution at supersonic 
speeds is presented in a double-integral form depending ex- 
plicitly upon the source-distribution function. The drag 
integral is applicable to a larger class of bodies than the usual 
slender-body drag expression and reduces to the Ward drag 
expression in the proper limit. Results for cones indicate 
that the drag integral with mass-flow continuity as the 
boundary condition gives good first-order results in the mod- 
erate supersonic range. 

The minimum-wave-drag problem for a transition section 
connecting two semi-infinite cylinders is solved with the aid 
of a suitably chosen auxiliary condition. The source distri- 
bution, the minimum drag, the slopes at the two ends of the 
section, and the radius at an intermediate point are obtained 
in terms of elementary functions. The entire shape is ob- 
tained in an integral form amenable to numerical evaluation. 
The minimum-drag shape obeys the Gothert similarity rule. 
The minimum- drag shape and the minimum drag are un- 
changed to the first order when the flow direction is reversed. 
The projectile tip is a special case and is compared with the 
Von KArmfin projectile tip. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 1, 1954- 



Figurb 7.- — Radius at intermediate point and slopes at two ends of 
minimum-extemal-drag transition section for five configurations. 
pRi— pRx~0.05. Entire contour shown for pRi + pih=0.118. 


APPENDIX 

INTEGRATION IN EQUATION (11) ° 


The integral of equation (11) is 

Z=j^ (s— g+a :—y)dx 




(Al) 


where the lower limit is the larger of the two quantities 
£4-/35! or !j4-/35i. The equation may be evaluated in the 
following manner. In order to convert the radicand in the 
denominator of the integrand into a quadratic function of the 
square of the integration variable, the variable is changed by 
using 

*=S4-^ (A2) 

Substitution of equation (A 2) into equation (Al) gives 

2 SdS 


/= 


s ° ^ s*—2 


(A3) 


where S L and So are the appropriate limits for the variable S. 
Since the numerator of the integrand in equation (A3) is an 
odd polynomial in S, considerable simplification is effected 
by taking a new variable z proportional to S 3 , such that 


s» =g+ (i_i) a +w 

Substituting equation (A4) into equation (A3) yields 


(A4) 


(AS) 


or 


K 


1= cosh -1 


(a ;—£)(x—ij)— 


WIT 


=cosh -1 


(t~v) 

(L-£)(£-n)-W 

05i(£— 4) 


IJt+wi 


(A6) 
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